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ABSTRACT 

Aims. We analyze the distribution of the interstellar matter in the environs of the Hn region Gum 31, excited by the open cluster 
NGC 3324, located in the complex Carina region, with the aim of investigating the action of the massive stars on the surrounding 
neutral material. 

Methods. We use neutral hydrogen 21cm-line data, radio continuum images at 0.843, 2.4 and 4.9 GHz, I2 CO(1-0) observations, and 
IRAS and MSX infrared data. 

Results. Adopting a distance of 3 kpc for the Hn region and the ionizing cluster, we have derived an electron density of 33±3 cirr 3 
and an ionized mass of (3.3±l.l)xl0 3 M Q based on the radio continuum data at 4.9 GHz. The Hi 21-cm line images revealed an Hi 
shell surrounding the Hn region. The Hi structure is 10.0± 1.7 pc in radius, has a neutral mass of 1500±500 M Q , and is expanding at 
11 kms -1 . The associated molecular gas amounts to (1.5±0.5)xlO s M Q , being its volume density of about 500 cirr 3 . This molecular 
material probably represents the remains of the cloud where the young open cluster NGC 3324 was born. The difference between the 
ambient density and the electron density of the H n region suggests that the H n region is expanding. 

The distributions of the ionized and molecular material, along with that of the emission in the MSX band A suggest that a photodis- 
sociation region has developed at the interface between the ionized and molecular gas. The copious UV photon flux from the early 
type stars in NGC 3324 keeps the H n region ionized. We conclude that either the massive stars in the open cluster have weak stellar 
winds or the stellar winds have blown during a very short period of time to create an interstellar bubble in an interstellar medium as 
dense as observed. 

The characteristics of a relatively large number of the IRAS, MSX, and 2MASS point sources projected onto the molecular envelope 
are compatible with protostellar candidates, showing the presence of active star forming regions. Very probably, the expansion 
of the H ii region has triggered stellar formation in the molecular shell. 

Key words. ISM: bubbles - H n regions - ISM: individual objects: Gum 3 1 - stars: early-type - stars: individual: HD 92206 



^ ■ 1. Introduction (PDR). Deharveng et al. (2005) have detected these PDRs in a 

„, . ,. r number of H ii regions. 

The interstellar medium associated with the birth place or mas- 

;3 sive stars, like O or early B-type stars, is made up of dense gi- In the present study, we analyze the distribution of the ion- 
ant molecular clouds. Massive stars are characterized by intense ized and neutral material associated with the H n region Gum 3 1 
photon fluxes and powerfull winds, which interact with their lo- (Gum 1955) based on H 1 21 cm line emission data, radio contin- 
cal medium ionizing and pushing the surrounding material. An uum information at different frequencies, and IR and molecular 
Hn region is a direct consequence of the high rate of Lyman data. 

continuum luminosity. At first time, the H n region is a small and ^ H n region Gum 3 { is about ,5, m size and approx . 

high density region, commonly named ultra-compact H 11 region imately circular m shape (Figure 1} It is i ocated at (l,b) = 

(UC H 11 Wood & Churchwell 1989). If the photon flux rate of (286 °12', _0°12') in the complex region of Carina and is consid- 

the massive star is sufficiently high, the H 11 region evolves into ered a member of the Car B 1 association. The SuperCOSMOS 

a noimal Hn legion. image (Parker et al. 2005) shows a quite inhomogeneous Hn 

A neutral shell encircles the H 11 region during the expanding region> with a sharp and bright rim towardfj lower galactic lon . 

phase (e.g. Spitzer 1978). The presence of these neutral shells gimdes and lower galactic latitudeSj looking fainter and more 

is observed in the Hi 21 cm line emission distribution (e.g. diffuse towards higher galactic longitudes. 
Deharveng et al. 2003). Molecular line studies have allowed the 

identification of molecular gas following the outer borders of Based on data of the radio recombination lines (RRL) H 1 09ff 
H 11 regions, indicating the presence of photodissociation regions a nd H 1 1 0a at 5 GHz, Caswell & Haynes (1987) found that the 
LSR velocity of the ionized gas in Gum 31 is -18 km s , sim- 



Send offprint requests to: C. Cappa ilar to the velocities of other Hn regions in the area of the Car 

* Member of Carrera del Investigador, CONICET, Argentina OBI association (Georgelin et al. 1986), and derived an electron 

** Prof. Virpi Niemela passed away on December 18th, 2006. temperature T e = 7100 K. They estimated a flux density S sghz 
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Table 1. Radio data: relevant parameters. 
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Radio continuum at 4.85 GHz 


Angular resolution 


5:0 


rms noise 


10 mjy beam -1 


Radio continuum at 2.4 GHz 


Angular resolution 


10:4 


rms noise 


12 mjy beam -1 


Radio continuum at 0.843 GHz 


Angular resolution 


43" x 51" 


rms noise 


1 mjy beam -1 


Hi data 


Synthesized beam 


2:4x2:1 


Number of channels 


256 


Velocity coverage 


(-190,+230) km s" 1 


Velocity resolution 


1.64 km s -1 


RMS noise level 


1.6 K 


12 CO(1-0) data 


Angular resolution 


2:7 


Velocity coverage 


(-50,+50) km s -1 


Velocity resolution 


0.2 km s -1 


RMS noise level 


1.0 K 



Fig. 1. SuperCOSMOS image showing the Hn region Gum 31. 
The cross marks the position of the multiple system HD 92206. 
The intensity units are arbitrary. 

The excitation sources of the Hn region Gum 31 are the 
OB star members of the open cluster NGC 3324. The bright- 
est star in this cluster is HD 92206, which is the visual double 
star IDS 10336-5806 in the Index catalogue (Jeffers et al. 1963) 
with a 1 mag fainter companion (HD 92206B) placed 5" to the 
East. Another bright cluster member is located 35" to the SW. 
This star is CD-57°3378, also referred to as HD 92206C in the 
literature. 

The three brightest stars in NGC 3324 have published spec- 
tral types. Both HD 92206A and B are classified as 06.5V, and 
HD 92206C as 09.5V by Mathys (1988). Walborn (1982) classi- 
fies HD 92206A as 06.5V(n) and the component C as 08.5Vp. 

Moffat & Vogt (1975) first carried out photometric observa- 
tions of about 12 stars in the cluster and estimated a color ex- 
cess E(B — V) ~ 0.45+0.05 mag and a distance d = 3.3 kpc. 
Claria (1977), using UBV photometry, confirmed previous re- 
sults by Moffat & Vogt and estimated d — 3.1 kpc. Vazquez & 
Feinstein (1990) derived a distance d = 3.6 kpc for the cluster 
from UBVRI photometry. More recently, Carraro et al. (2001) 
found about 25 new possible cluster members and derived d — 
3.0+0.1 kpc. On the other hand, distance estimates for Car OBI 
are in the range 1.8-2.8 kpc (Walborn 1995). Bearing in mind 
these results we adopted d = 3.0+0.5 kpc for both Gum 31 and 
the ionizing cluter. 

Carraro et al. (2001) find evidence for pre-main sequence 
members beginning at about late B spectral type, which suggests 
an extremely young age for NGC 3324 (< (2-3)xl0 6 yr). The O- 
type members would be stars recently arrived on the Zero Age 
Main Sequence. 

2. Data sets 

2.1. Radio data sets 

We have analyzed the radio continuum emission in the region 
of Gum 31 using data obtained at 0.843, 2.42 and 4.85 GHz, 
which were extracted from the Sydney University Molonglo Sky 



Survey (SUMSS) (see Sadler & Hunstead 2001 for a descrip- 
tion), the survey by Duncan et al. (1995), and the Parkes-MIT- 
NRAO (PMN) Southern Radio Survey (see Condon et al. 1993 
for a complete description of this survey), respectively. 

We have used Hi data from the Southern Galactic Plane 
Survey (SGPS) to analyze the neutral gas distribution in the 
environs of Gum 31. These data have been obtained with the 
Australia Telescope Compact Array (ATCA) and the Parkes 
Radiotelescope (short spacing information). A Hanning smooth- 
ing (Rohlfs 1986) was applied to the Hi data to improve the 
signal to noise ratio. A description of this survey can be found in 
McClure-Griffiths et al. (2005). 

The distribution of the molecular material in the region was 
studied using 12 CO(1-0) line data at 115 GHz obtained with the 
NANTEN 4 m telescope of Nagoya University at Las Campanas 
Observatory of the Carnegie Institution of Washington, and pub- 
lished by Yonekura et al. (2005). 

The main observational parameters of these data bases are 
listed in Table 1 . 

2.2. Infrared data 

We have also investigated the dust distribution using high- 
resolution (HIRES) IRAS, and MSX data obtained through 
/PACfl The IR data in the IRAS bands at 60 and 100/mi have 
angular resolutions of LI and Y.9. The images in the four MSX 
bands (8.3, 12.1, 14.7, and 21.3 /mi) have an angular resolution 
of 1873. MSX units were converted into MJy ster -1 by multiply- 
ing the original figures by 7.133xl0 6 (Egan et al. 1999). 
No Spitzer data exist on this region. 

3. Results 

3.1. Radio continuum emission 

Figure 2 displays the radio continuum image at 843 MHz. The 
image reveals a radio source of 15' in size, coincident in position 
with the H 11 region. There is a remarkable correlation between 

1 IPAC is funded by NASA as part of the IRAS extended mission un- 
der contract to Jet Propulsion Laboratory (JPL) and California Institute 
of Technology (Caltech). 
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Fig. 2. Top: Radio continuum image at 843 MHz. The grayscale 
is from -5 to 180 mJy beam -1 . The contour lines are 5 (5<x), 
20, 50, 100, 150, and 200 mJy beam -1 . The cross indicates 
the position of HD 92206. Bottom: Overlay of the image at 843 
MHz (contour lines) and the superCOSMOS image of Gum 3 1 
(grayscale). 



radio and optical emission regions. The strongest radio emission 
region coincides with the brightest section of the optical rim at / 
=s 286° 10'. The region of diffuse emission near (l,b) = (286°20', 
-0°4') has also a weak radio counterpart. The fainter optical re- 
gions at (l,b) = (286°13', -0°15') and (286°20', -0°7') also cor- 
relate with weak radio emission. Along with the optical image, 
this radio image shows that the H n region is far from being ho- 
mogeneous and that it is non-uniform in density. Gum 3 1 is also 
detected in the surveys at 4.85 and 2.4 GHz as an isolated radio 
source. The flux densities at 4.85 and 2.42 GHz are 37.7+2.5 
Jy (coincident with the previous estimate by Caswell & Haynes 
[1987]) and 40.7+2.9 Jy, respectively. 

The spectral index of the radio source, as derived from the 
emissions at 2.4 and 4.85 GHz, is a = -0.09+0.20, confirming 
its thermal nature. This result is compatible with the detection of 
RRL at 5 GHz. 



3.2. Hi results 

Circular galactic rotation models predict negative radial veloc- 
ities of up to about -12 km s _1 in the line of sight towards I 
= 286°. However, radial velocities observed in this section of 
the Galaxy are more negative (of up to -30 km s _1 , Brand & 
Blitz 1993), indicating the presence of non-circular motions. 
Consequently, we paid special attention to the neutral gas emis- 
sion distribution at negative radial velocities of up to -50 km 
s- 1 . 

Figure 3 displays a series of H i images within the velocity 
interval from -12.8 to -32.5 km s _1 in steps of 3.3 km s _1 . The 
presence of a region of low H i emission surrounded by regions 
of enhanced emission approximately centered at the position of 
HD 92206 is clearly identified within the velocity range -12.8 to 
-32.5 km s _1 . An almost complete envelope encircles the void. 

The top panel of Figure 4 displays the H i emission distribu- 
tion within the velocity range -29.3 to -16.5 km s _1 , where the 
relatively thick envelope, of about 8'-10' (or 7.0-8.7 pc at d = 3 
kpc), is better defined. 

The bottom panel of Figure 4 displays an overlay of the H i 
emission distribution and the SuperCOSMOS image of Gum 3 1 . 
The H i envelope clearly anti-correlates with the ionized nebula. 
This envelope is less bright towards higher galactic latitudes. 

The systemic velocity of the structure, which corresponds to 
the velocity at which the feature presents its largest dimensions 
and deepest temperature gradient, is about -23 km s _1 . This ve- 
locity is similar to the velocity of the ionized gas (-18 km s _1 , 
see §1) as obtained from radio recombination lines, within the 
errors. 

The morphological agreement between the optical nebula 
and the H i shell and the agreement in velocity between the ion- 
ized and neutral materials indicate that the H i feature is the neu- 
tral counterpart of the ionized nebula. 

3.3. The CO emission distribution 

The distribution of the molecular gas is shown in Fig. 5. The 
upper panel displays the 12 CO(1-0) emission distribution in 
grayscale and contour lines, while the bottom panel displays 
an overlay of the CO contour lines and the optical image of 
Gum 3 1 . The CO gas distribution displayed in the figure was 
obtained from the molecular data cube kindly provided by Y. 
Yonekura. We integrated the 12 CO(1-0) emission within the 
velocity interval -27.2 to -14.0 km s This velocity range is 
slightly different from the one used by Yonekura et al. (2005) 
(-30 to -10 km s '), since no CO emission was detected for 
velocities v < -27 km s and v > -14 km s 1 associated with 
Gum 31. Both images, the one by Yonekura et al. and the one 
in Fig. 5 are essentially the same. 

Intense CO emission regions encircle the brightest sections 
of the optical nebula with CO clumps strikingly bordering the 
bright rim at (286° 10', -0°10'), and near (286°23', -0°15'), 
where the nebula appears diffuse. The CO envelope is open to- 
wards (l,b) = (286°25', -0°5'). The region of relatively faint op- 
tical emission near (l,b) = (286°20\ -0°4') coincides with a low 
CO emission region. 

The close morphological agreement between the optical and 
the molecular emissions in the brightest optical region indicates 
that the molecular material is interacting with Gum 3 1 . 

The thick lines in the bottom panel of Fig. 5 delineate the 
C 18 cores found by Yonekura et al. (2005). The dense cores 
coincide with the the brightest 12 CO emission regions. 
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Fig. 3. Series of Hi brightness temperature images for the velocity interval -12.8 to -32.5 km s _1 averaged over 3.3 km s _1 . The 
greyscale is from 70 to 120 K for the images at v > -26.0 km s _1 , and from 10 to 60 K for the images at v < -26.0 km s _1 . The 
contour lines are from 10 to 1 10 in steps of 10 K. The star indicates the position of HD 92206. The velocity interval corresponding 
to each image is indicated. 



The comparison of the molecular gas distribution as shown 
by the 12 CO emission (Fig. 5) with the Hi structure (Fig. 4) 
shows that the neutral envelope around the Hn region has a 
molecular component. Although, both the H i and the molecular 
envelopes are approximately coincident, the regions of brightest 
molecular emission anticorrelate with the H i maxima. 

The comparison between the ionized, neutral atomic, and 
molecular distributions around Gum 31 suggests the presence 
of a prominent photodissociation region (PDR) bordering the 
brightest ionized regions. Very probably, the CO emission shows 
the remaining parental molecular material. 



3.4. The emission in the infrared 

The distribution of the IR emission at 60 and 100 yum, due to ther- 
mal dust emission, is displayed in Figure 6. The upper panels 
show the emission at both wavelengths while the bottom pan- 
els display overlays of the IR emission and the optical image 
of Gum 3 1 . The images reveal an IR structure which is brighter 
near (l,b) = (286°23',-0° 15') and where the optical nebula has its 
sharpest border, and weaker on the opposite side. The brightest 
IR emission regions at 60/im and lOO/im are projected onto the 
neutral envelope, delineating the ionized nebula. The IR emis- 
sion at (286°23',-0°15') detected at both wavelengths coin- 



cides with strong 12 CO(1-0) emission. 12 CO emission also ap- 
pears bordering the two IR clumps located near the bright 
rim at (l,b) = (286°10', -0°10) detected at 60^m, and in be- 
tween. The distribution of the IR emission at both wavelengths 
shows the presence of dust most probably related to the sur- 
rounding Hi and molecular shells. 

Following the procedure described by Cichowolski et al. 
(2001), we derived the color temperature of the dust associated 
with the H n region and the neutral envelope based on the IR 
fluxes at 60yum and 100/im. Taking into account different val- 
ues for the background emission, we found Tj = 34+7 K. The 
range of temperatures corresponds to n = 1-2 and to different 
IR background emissions. The parameter n is related to the dust 
absorption efficiency (k v oc y". We adopted k v = 40 cm 2 g 
This value was derived from the expressions by Hildebrand 
(1983) for n =1-2 in the range 60 and lOO^rn. The obtained 
dust temperature is typical for H n regions. 

Figure 7 shows an overlay of the distribution of the emission 
in the MSX bands A and E, and the optical image. The emission 
in band A closely follows the brightest sections of the nebula, 
and correlates with the neutral gas. Particularly, the brightest 
regions emitting at 8.3 /jm coincide with the the dense cores 
2 and 6 (see Fig. 5) found by Yonekura et al. (2005). On the 
contrary, the emission in band E appears clearly associated with 
the ionized gas. Note that the strongest emission region in band 
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Fig. 6. Top: Far infrared IRAS images at 60 and 100//m towards Gum31. The grayscale is from 300 to 3000 MJy ster~' for the 
image at 60/im, and from 600 to 4500 MJy ster -1 for the image at 100//m. The contour lines are from 200 to 1000 MJy ster~' in 
steps of 200 MJy ster 1 and from 1000 to 3000 MJy ster 1 in steps of 500 MJy ster 1 for the image at 60//m, and from 400 to 1000 
MJy ster~' in steps of 200 MJy ster~' and from 1000 to 3000 MJy ster -1 in steps of 500 MJy ster" 1 for the image at lOOyum. Bottom: 
Overlay of the optical (SuperCOSMOS) (grayscale) and the IR images. 
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Fig. 7. Overlay of the MSX infrared images (contours) corresponding to bands A (8.28/im) and E (21.34 /urn), and the 
SuperCOSMOS image of the nebula (grayscale). The contour lines are 25, 39, 57, 85, 114, and 140 MJy ster 1 for band A; 
and 36, 46, 57, 85, 1 14, and 140 MJy ster 1 for band E. 
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Fig. 9. Top: Point sources from the IRAS (triangles), MSX (circles), and 2MASS catalogues overlaid onto the image shown in the 
bottom panel of Fig. 5. 



A at (l,b) = (286°23',-0°15') coincides with bright molecular 
and far infrared clumps. 

The emission distribution in band A is most probably related 
to emission from polycyclic aromatic hydrocarbons (PAHs). 
According to Cesarsky et al. (1996), these dust grains can not 
survive inside the Hn region, but on the neutral PDR, where 
they radiate in the PAH bands at 7.7 and 8.6 /im, included in the 
MSX band A. MSX band E, on the contrary, includes contin- 
uum emission from very small grains, which can survive inside 
ionized regions (cf. Deharveng et al. 2005), and an important 
contribution from nebular emission lines. 

In summary, the distribution of the neutral atomic and molec- 
ular gas, and that of the interstellar dust reveals the presence of 
a neutral shell surrounding the Hn region. The distribution of 
the molecular gas and that of the emission in the MSX band A 
related to PAHs strongly suggests the presence of a PDR at the 
interface between the ionized and the molecular gas. 

Both the infrared emission at 60 and lOOyum and the optical 
image suggest that the faint optical emission region near at (l,b) 
= (286°20', -0°4') (indicated by an arrow in the bottom right 
image of Fig. 6) are also linked to Gum 3 1 . As pointed out in 
Sect. 3.1, this region is also faint in the radio continuum (see Fig. 
2), and corresponds to weak regions in the H i and CO envelopes 
suggesting that the ionizing flux and the stellar wind energy of 
the massive stars in the open cluster may drain through these 
regions to the general ISM. The situation resembles the case of 
the stellar wind bubble around WR 23 (Cappa et al. 2005). 



4. Discussion 

4.1. Physical parameters 

The main parameters of the dust and the ionized and neutral gas 
related to Gum 3 1 are summarized in Table 2. The parameters of 
the ionized gas were derived from the image at 4.85 GHz. The 
uncertainty in the flux density corres ponds to an error of 0.1 
Jy beam -1 in the estimate of the radio continuum background. 
The electron density and the H n mass were obtained from the 
expresions by Mezger & Henderson (1967) for a spherical Hn 
region of constant electron density (rms electron density n e ). The 
presence of He n was considered by multiplying the H i mass by 
1.27. The number of UV photons necessary to ionize the gas 
log Niy-c was derived from the radio continuum results. Errors 
in the linear radius, in the rms electron density, and in the exci- 
tation parameter come from the distance uncertainty. The high 
electron density of the H n region is compatible with the rela- 
tively short lifetime for the cluster. 

The parameters of the neutral atomic gas includes: the (l,b) 
position of the centroid of the Hi shell, the velocity interval 
spanned by the structure, the systemic and expansion velocities, 
the radius of the neutral gas structure, and the associated atomic 
mass. 

The expansion velocity was estimated as in previous papers 
(see Cappa et al. 2005 and references therein) as v exp = (v2-vi)/2 
+ 1 .6 km s _1 . The extra 1 .6 km s _1 allows for the presence of H i 
in the caps of the expanding shell, which are not detected in the 
present case. 
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Adopted distance (kpc) 


3.0+0.5 


Hn region 




Flux density at 4.85 GHz (Jy) 


37.7±2.5 


Angular radius (arcmin) 


7.5+0.2 


Linear radius (pc) 


6.5±1.0 


Emission measure (pc cirr 6 ) 


(1.5±0.2)xl0 4 


rms electron density [« e ](cirr 3 ) 


33±3 


Used Lyman UV photons [logN Ly - c ] (s _1 ) 


49.0±0.1 


Ionized mass (M Q ) 


3300±1100 


Neutral atomic shell 




(l,b) centroid 


286° 15' -0° 10' 


Velocity interval Vi, V2 (km s -1 ) 


-13,-32 


Hi systemic velocity (km s -1 ) 


-23 


Expansion velocity v exp (km s -1 ) 


11 


Radius of the H i structure (arcmin) 


11.5 


Radius of the H i structure R (pc) 


10.0±1.7 


Atomic mass in the shell (M Q ) 


1500+500 


Molecular shell 




Velocity interval Vi,V2 (km s -1 ) 


-27.2,-14.0 


Mean H2 column density (cm -2 ) 


1.2xl0 22 


Molecular mass (M Q ) 


(l.l±0.5)xl0 5 


Dust related to the Hn region and the neutral shell 

Total dust mass (M ) 60+20 

Dust color temperature (K) 34+7 



The radius of the Hi structure was estimated from Fig. 4 
and corresponds to the position of the maxima in the shell. The 
neutral atomic mass corresponds to the mass excess in the shell, 
assuming that the gas is optically thin and including a He abun- 
dance of 10%. 

The H2 column density (Nm) and the molecular mass were 
estimated from the 12 CO data, making use of the empirical rela- 
tion between the integrated emission Wco (= jTdv) and Nhi- 
We adopted N H2 = (1.06±0.14) xW C0 x 10 20 cirr 2 (K km 
s -1 ) -1 , obtained from y-ray studies of molecular clouds in the 
Orion region (Digel et al. 1995). To derive the molecular mass, 
we integrated the CO emission within a circle of about 16' in 
radius (= 15 pc at 3.0 kpc), centered at (l,b) = (286°16',-0°8'). 

The ambient density derived by distributing the molecular 
mass over a sphere of 15.0 pc in radius is ^500 crrT 3 , reinforcing 
the idea that the molecular gas represents the remains of the orig- 
inal material where the open cluster NGC 3324 was born. The 
difference between the electron density and the ambient density 
indicates that the H 11 region is expanding. 

The optical appearance of the H 11 region, almost spherical 
without clear evidences of a central cavity, suggests that the mas- 
sive stars in the cluster have weak stellar winds (which is sup- 
ported by the short age of the cluster) and/or have existed during 
a very short period of time to create an interstellar bubble in an 
interstellar medium as dense as observed here. 

The sketch of Fig. 8 shows schematically the distribution 
of the different gas components. Although Hi and molecular 
clumps anticorrelate in position, the large scale H 1 and molecu- 
lar gas coincide in position all around the nebula but in the region 
at higher galactic latitudes, where the CO emission is observed 
outside the H 1 region. The large amount of molecular gas com- 
pared to the neutral atomic gas supports the idea that the CO 
emission that encompassed the ionized nebula represents the re- 
mains of the molecular material where the open cluster was born. 
The fact that part of this H 1 gas may have been originated in the 
photodissociation of the molecular gas. 



4.2. Energy budget 

Bearing in mind that massive stars have a copious UV flux ca- 
pable of ionizing the surrounding H 1 gas we investigate in this 
section whether the massive stars in NGC 3324 can provide the 
energy to ionize the gas. 

As described in Sect. 1, the only O-type stars in the open 
cluster are the multiple system HD 92206, which contains three 
O-type stars classified as 06.5V, 06.5V and 08.5 (Walborn 
1982, Mathys 1988). Considering that HD92206A is almost 1 
mag brighter than component B (Claria 1977), the similar spec- 
tral types found by Mathys (1988) suggest that HD92206A is 
probably a spectroscopic binary 06.5V + 06.5V. Taking into 
account the UV photon fluxes emitted by the stars, N, (s -1 ), de- 
rived by Martins et al. (2005) from stellar atmosphere models, 
a group of four massive stars having the spectral types indicated 
above have a total UV photon flux corresponding to log N* = 
49.4. By comparing this value with the UV photons used to ion- 
ized the gas listed in Table 2, we conclude that the massive O- 
type stars in NGC 3324 can maintain the H 11 region ionized. 

The radius of the Stromgren sphere formed in a region with 
the volumetric ambient density listed in Table 2 is of ^ 5.0 pc, 
lower than the radius of the ionized region. Following Spitzer 
(1978), we estimated that the Hn region have been expanding 
during as 2xl0 5 yr. The ambient density we have derived from 
the 12 CO(1-0) data is larger than the rms electron density, 33 
cirr 3 (see Table 2), giving additional support for the interpreta- 
tion that the H 11 region is expanding. 



5. Stellar formation 

Stellar formation may be induced in the surrounding molecu- 
lar envelope, where the presence of high density regions favors 
the conditions which lead to form new stars. Among the differ- 
ent processes that induce star formation, the collect and collapse 
process proposed by Elmegreen & Lada (1977) may work effi- 
ciently in the dense molecular shells around H 11 regions. 
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Table 3. IR Point sources and YSOs candidates from the IRAS, MSX and 2MASS catalogues 



# 


l[0 '] 


b[o '] 


IRAS source 




Fluxes [Jy] 




L IRAS [WL Q ] 










12/jm 


25//m 


60//m 


100/jm 




1 


285°59:82 


+0°05:82 


10349-5801 


1.6 


2.5 


27.9 


104 


1.5 


2 


286°04:98 


-0°25:08 


10335-5830 


0.8 


3.5 


22.8 


111 


1.5 


3 


286°08:58 


-0°06:96 


10351-5816 


6.3 


6.8 


315 


1480 


18 


4 


286°08:64 


-0° 18:78 


10343-5826 


5.5 


4.3 


103 


340 


5 


5 


286° 10:92 


-0°14:7 


10349-5824 


5.9 


9.0 


163 


1660 


18 


6 


286 12! 18 


+0 10' 20 


10365-5803 


7.2 


86 


1170 


2780 


43 


/ 


ZOO lj.12 






t i 
1.1 


3.3 


11.7 


1430 


1 A 

14 


8 


286 17:34 


+0 00:42 


10365-5814 


2.3 


1.7 


85 


246 


4 


9 


286 22:5 


-0°15:3 


10361-5830 


12.4 


38.4 


626 


2160 


30 


10 


286 26!94 


~0°22:92 


10361-5839 


2.8 


5.4 


51.6 


307 


4 


11 


286 29:34 


-0° 16:68 


10368-5835 


2.1 


2.3 


84 


271 


4 


12 


286 33:72 


ao Any ad 

-0 07:08 


10379-5828 


1.0 


2.0 


13.8 


540 


5 


# 


/[C] 


b[o '] 


MSX source 




Fluxes [ J Y] 














8/jm 


12//m 


14yum 


21^111 


Class. 


13 


286°09:78 


-o° 11:28 


G286. 1626-00. 1877 


0.7311 


1.298 


1.256 


2.806 


CHu 


14 


286° 12!36 


-0°09:66 


G286.2056-00.1611 


0.1585 


0.9313 


2.151 


6.904 


MYSO 


i c 


2S6 12.45 


no 1 n' qo 
— U 10.32 


O280.2U / /-U0. 1 /2U 


U.OS /U 


0.9481 


1.764 


2.y 


CM II 


16 


286° 12:54 


+o°io:i4 


G286. 2086+00. 1 694 


1.353 


2.882 


7.182 


40.57 


MYSO 


17 


286° 12:6 


-0° 10:68 


G286. 2096-00. 1775 


0.2202 


0.9538 


1.368 


7.206 


CHu 


18 


286° 15:42 


-0° 19:44 


G286. 2566-00.3236 


2.04 


2.151 


1.34 


4.328 


CHu 


19 


286°2i:48 


-0° 17:58 


G286. 3579-00.2933 


0.7126 


1.815 


2.677 


6.065 


MYSO 


20 


286 22:5 


-0° 15:78 


G286. 3747-00.2630 


3.591 


4.756 


2.409 


7.577 


CHu 


21 


286°22:62 


-0° 15:36 


G286. 3773-00.2563 


1.628 


2.918 


3.855 


12.1 


CHu 


# 


l[o'] 


b[o>] 


2MASS source 


-/[mag] 


#[mag] 


A:, [mag] 


(J-H) 


(H-K) 


22 


285°54:42 


-o° 13:74 


10350210-5831039 


10.584 


10.631 


10.564 


-0.047 


0.067 


23 


286°04:44 


-0°00:24 


10365972-5824186 


12.042 


10.51 


9.284 


1.532 


1.226 


24 


286°06:36 


-0°08:64 


10364112-5832326 


11.367 


11.385 


11.28 


-0.018 


0.105 


25 


286°06:54 


-0°08:40 


10364296-5832267 


10.465 


10.513 


10.439 


-0.048 


0.074 


ZD 


Zoo Uo.Zz 


n o O/i ' ^A 
— U U4.JO 


lUi /Usyo-jB2yjj J 


12. / /D 


11.622 


10.896 


i i 

i.ijj 


n toa 
U. /20 


27 


286°09:60 


-0=00:30 


10373406-5826540 


12.487 


11.454 


10.513 


1.033 


0.941 


28 


286°09:78 


-o° 11:22 


10365396-5836293 


12.242 


10.99 


10.103 


1.252 


0.887 


29 


286°10:08 


-0°09:i2 


10370395-5834489 


10.772 


9.712 


8.685 


1.06 


1.027 


30 


286° 10:20 


-o°n:io 


10365749-5836366 


13.99 


12.815 


11.984 


1.175 


0.831 


31 


286°li:i6 


-0°02:64 


10373574-5829405 


15.318 


13.343 


11.614 


1.975 


1.729 


32 


286° 12: 18 


-0°02:i6 


10374424-5829451 


13.004 


11.588 


10.685 


1.416 


0.903 


33 


286° 12:96 


-0°04:86 


10373956-5832311 


10.926 


10.752 


10.554 


0.174 


0.198 



To investigate the presence of protostellar candidates in this 
region we used data from the IRAS, MSX and 2MASS point 
source catalogues. We searched for point sources in a region of 
about 20' in radius centered at the position of the NGC 3324. 



5.1. IRAS sources 



The IRAS point source catalogue allows identification of pro- 
tostellar candidates following the criteria by Junkes et al. 
(1992). Sources with quality factors Q h o + 2ioo > 4 were 
considered. Twelve out of the thirteen point sources have IR 
spectra compatible with protostellar objects. The names of 
the IRAS protostellar candidates, their (l,b) position, their 
fluxes at different IR wavelengths, and their luminosity de- 
rived following Yamaguchi et al. (2001), along with a ref- 
erence number are listed in Table 3. The IRAS protostellar 
candidates are indicated as filled triangles in Fig. 9, superim- 
posed onto the molecular gas and ionized gas distributions. 
Each source can be identified in the figure by its reference 
number. They are probably Class objects. 



5.2. MSX sources 

Massive young stellar objects (MYSOs) can be identified 
from the MSX point catalogue following the criteria by 
Lumsden et al. (2002). A total of 310 MSX point sources were 
found to be projected onto the area. 

Lumsden et al. (2002) found that MYSOs have infrared 
fluxes with ratio F 2 i/F% > 2 and F u /Fu > 1, where F 8 , F\ 2 , 
F 14, and F 2i are the fluxes at 8.3, 12, 14, and 21 /on. For 
compact Hn regions, ratios are F 21 /F% > 2 and Fu/Fn < 1. 
Taking into account sources with flux quality q > 2, we were 
left with 14 sources after applying Lumsden et al.'s criteria 
indicated above. 

Three out of the 14 sources can be classified as MYSOs, 
while 6 sources are compact H n regions, also indicative of 
active stellar formation. The nine sources are listed in Table 
3 and indicated in Fig. 9 as circles. 

5.3. 2MASS sources 

Point sources with infrared excess, which are candidates to 
young stellar objects, were searched for in the 2MASS cata- 
logue (Cutri et al. 2003), which provides detections in three 
near IR bands: J, H, and K s , at 1.25, 1.65, and 2.17yum, 
respectively. A total of 2x 10 4 sources are projected onto 



Table 3. Cont. 
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# 


H°'] 


b[o>] 


2MASS source 


-/[mag] 


W[mag] 


K s [mag] 


(J-H) 


(H-K) 


34 


286° 13:02 


-0° 10:86 


10371717-5837460 


11.793 


11.773 


11.665 


0.02 


0.108 


35 


286°13:32 


-0°08:46 


10372824-5835492 


12.191 


11.301 


10.448 


0.89 


0.853 


36 


286° 1338 


-o° 10:20 


10372226-5837229 


7.563 


7.588 


7.479 


-0.025 


0.109 


37 


286° 13:92 


-0° 17:87 


10365763-5844052 


12.22 


11.882 


11.51 


0.338 


0.372 


38 


286° 15:54 


-0° 19:44 


10370125-5846295 


11.446 


11.03 


10.72 


0.416 


0.31 


39 


286° 16:92 


-0°0i:38 


10381935-5831264 


12.613 


11.905 


11.42 


0.708 


0.485 


40 


286° 19:26 


-0°06:84 


10381421-5837192 


12.635 


11.886 


11.349 


0.749 


0.537 


41 


286° 19:62 


-0°06:84 


10381639-5837318 


12.71 


11.744 


11.11 


0.966 


0.634 


42 


286°20:i0 


-0° 19:80 


10373105-5849026 


15.155 


12.8 


11.127 


2.355 


1.673 


43 


286°2i:06 


-0° 16:86 


10375219-5847133 


12.271 


11.363 


10.675 


0.908 


0.688 


44 


286°2i:90 


-0° 04:92 


10383875-5836566 


12.116 


11.636 


11.187 


0.48 


0.449 


45 


286°22'92 


-0° 13:20 


10381461-5844416 

1 v JU l^V 1 -J O III 1 \J 


12.367 


11.815 


11.398 


0.552 


0.417 


46 


286° 24: 92 


-0° 18:66 


10380736-5850240 


12.477 


11.538 


10.78 


0.939 


0.758 


47 


286°28:il 


-0°24:50 


10380702-5857039 


12.724 


11.599 


10.726 


1.125 


0.873 


48 


286°28:44 


-0°23:71 


10381226-5856318 


13.348 


12.274 


11.355 


1.074 


0.919 


49 


286°29:46 


-0°06:30 


10392451-5841486 


13.836 


12.561 


11.619 


1.275 


0.942 


50 


286°29:46 


-o°08:io 


10391799-5843257 


13.195 


12.322 


11.624 


0.873 


0.698 


51 


286°30:90 


-0°06:30 


10393410-5842321 


11.701 


10.423 


9.334 


1.278 


1.089 


52 


286° 3 i:27 


-0°28:40 


10381363-5902003 


13.368 


12.557 


11.982 


0.811 


0.575 



a circular region of 20' in radius. We took into account 
sources with S /N > 10 (corresponding to quality "AAA"). 
Only sources with K s < 12 were included. The last criterium 
corresponds to stars with spectral types earlier than B3 at 
a distance of 3.0 kpc. Following Comeron et al. (2005) and 
Romero (2006), we determined the parameter q = (J - H) - 
1 .83 x (H-K s ). Sources with q < -0.15 are classified as objects 
with infrared excess that may be young stellar objects. After 
applying this criterium we were left with 31 sources, which 
are shown in the magnitud-color diagram of Fig. 10. The lo- 
cation of these sources is also marked in Fig. 9 as crosses. The 
main data of these sources are shown in Table 3. 

Although the diagram is not conclusive in identifying 
young stellar objects, the strong infrared excess of the 
sources is compatible with protostellar candidates. 

5.4. Distribution and characteristics of the YSOs 

Figure 9 shows that the IRAS and MSX point sources, and 
most of the 2MASS point sources classified as YSOs appear 
bordering the ionized region, projected onto the molecular 
envelope detected in 12 CO emission, close to the periphery 
of the H 11 region. Some of them are also coincident with the 
dense cores found by Yonekura et al. (2005) in C ls O. 

We will analyze some particular regions. 

The IRAS source #9, the MSX sources #19, #20, and #21, 
and the 2MASS sources #43 and #45, are projected onto a 
12 CO clump and onto the dense core 6 found in C ls O. The 
presence of these sources indicates that stellar formation is 
going on in this particular molecular clump. As suggested 
by Yonekura et al., a star cluster including massive stars is 
probably being formed in this region. 

Particularly interesting is the bright rim region at (l,b) 
* (286°10',-0 o 10). The 2MASS sources #28, #29, and #30, 
and the MSX source #13 are located in this area, projected 
onto the molecular envelope. IRAS sources #3 and #5 are also 
placed close to the border of the H 11 region, coincident with 
molecular material. Source #3 coincide with a region emit- 
ting in the MSX band A and with core 2. 

A bunch of protostellar objects is almost coincident in 
position with the open cluster NGC 3324: one MYSOs (MSX 



source #14), two compact H 11 regions (MSX sources #15 and 
#17), and the 2MASS sources #34 and #36. Some of these 
sources coincide with the loose IR cluster IC 2599 listed by 
Dutra et al. (2003). These facts suggest that stellar forma- 
tion is still going on in the region of NGC 3324, as previously 
found by Carraro et al. (2001). 

Also IRAS sources #6 and #10 are projected onto the 
dense cores 5 and 7, respectively. 

The large luminosities L mAS estimated for some of these 
sources suggest that they are protostellar candidates for mas- 
sive stars or star clusters. Also note that IRAS protostellar 
candidates projected onto the molecular envelope close to the 
periphery of the H 11 region have large luminosities in com- 
parison with most of the sources located far from the ion- 
ization front. This is compatible with previous findings by 
Dobashi et al. (2001). 

The sources IRAS 10355-5828 (286°16:92,-0°15:6) and 
MSX G286.2868-00.2604 (neither included in Table 3 nor 
shown in Fig. 9) coincide with HD 92207, a red supergiant 
whose membership to the open cluster is a matter of de- 
bate (Carraro et al. 2001). Its MSX fluxes correspond to an 
evolved object. 

The detection of protostellar candidates in the IRAS, 
MSX, and 2MASS data bases, strongly indicates that active 
stellar formation is currently going on in the molecular shell 
around Gum 31. 

To summ up, most of the protostellar candidates detected 
towards Gum 31 appear projected onto the shell detected in 
the 12 CO line, and coincide with the dense cores detected in 
C 18 emission. Some of them are located close to the bright 
sharp borders of the H 11 region, and near the open cluster 
NGC 3324. 

The presence of protostellar objects onto the molecular 
envelope bordering the ionized region indicates that star for- 
mation has been triggered by the expansion of the H 11 region. 
The distribution of the molecular and Hi gas around the ion- 
ized region suggests that star formation could be due to the 
collect and collapse process. 
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Fig. 4. Top: H i brightness temperature image corresponding to 
the velocity interval -29.3 to -16.5 km s _1 showing the H i struc- 
ture related to Gum 31. The grayscale is from 40 to 100 K. The 
contour lines are from 40 to 100 K in steps of 10 K. Bottom: 
Overlay of the SuperCOSMOS (grayscale) and Hi (contours) 
images showing the close correspondence of neutral and ionized 
gas emissions. 



6. Summary 

We have analyzed the interstellar medium in the environs of the 
H ii region Gum 3 1 to investigate the action of the massive stars 
in the exciting open cluster NGC 3324 on the surrounding neu- 
tral material. 

We based our study on H 1 21cm line emission data belonging 
to the Southern Galactic Plane Survey, radio continuum data at 
0.843, 2.4 and 4.85 GHz from the PMN Southern Radio Survey, 
12 CO data from Yonekura et al. (2005), and IRAS (HIRES) and 
MSX data. 

Adopting a distance of 3.0+0.5 kpc, we have derived an ion- 
ized mass of 3300±1100 M Q and an electron density of 33+3 
cirT 3 . The four O-type stars in the HD 92206 multiple system 
can provide the necessary UV photon flux to maintain the H n 
region ionized. 

The Hi emission distribution in the environs of Gum 31 
shows the presence of an H i shell approximately centered at the 
position of the multiple system HD 92206. The H i shell closely 
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Fig. 5. Top: 12 CO(1-0) emission distribution corresponding to 
the velocity interval -27.2 to -14.0 km s" 1 . The grayscale 
is from 15 to 150 K km s _1 . The contour lines are 20.2, 
40.5, 60.7, 81.0 and 101.2 K km s -1 . Bottom: Overlay of 
the SuperCOSMOS (grayscale) and the CO (contours) im- 
ages. The thick lines delineate the C ls O cores described by 
Yonekura et al. (2005). The crosses mark the core positions 
indicated in table 2 by Yonekura et al . 



encircles the optical nebula. It is detected within the velocity 
range -32 to -13 km s _1 and its systemic velocity of -23 km s _1 
is coincident, within errors, with the velocity of the ionized gas 
in the nebula (-18 km s _1 ). The Hi structure is 10.0+1.7 pc in 
radius and expands at about 11 km s _1 . The associated atomic 
mass is 1500+500 M . 

Molecular gas with velocities in the range -27.2 and -14.0 
km s~' surrounds the brightest parts of Gum 31. The sharp 
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CO clump 



Fig. 8. Schematic view of Gum 3 1 and the respective locations 
of the different gas components and of interstellar dust. 
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Fig. 10. Magnitud-color diagram showing 2MASS sources 
whith infrared excess. 



interface between the ionized and molecular material indicates 
that these gas components are interacting. We have estimated 
a molecular mass of (l.l+0.5)xl0 5 M Q . We suggests that the 
molecular gas represents the remains of the natal cloud 
where NGC 3324 formed. 

The volume density of the molecular cloud, higher than the 
electron density (33 cm 3 ), implies that the H n region is expand- 
ing. 

The emission in the far infrared correlates with the H n re- 
gion and the Hi envelope, indicating that the observed emis- 
sion is probably related the to the neutral and molecular en- 
velopes, and molecular envelopes. 

The distribution of the emission in the MSX band A, that 
closely delineates the H n region and correlates with the molec- 
ular emission, suggests the presence of a PDR at the interface 
between the ionized and molecular gas. 



A number of MSX, IRAS, and 2MASS point sources with 
IR spectra compatible with protostellar objects appear projected 
onto the molecular envelope, implying that stellar formation is 
active in the higher density cores of the molecular envelope 
around Gum 3 1 , where massive stars or star clusters are prob- 
ably being formed. 

The optical image of the nebula does not show clear evi- 
dences of a central cavity, as expected in a stellar wind bubble, 
suggesting that the massive stars in the cluster have weak stellar 
winds or have existed during a short period of time to develop 
an interstellar bubble in a high density interstellar medium. 

Acknowledgements. Prof. Virpi Niemela passed away a few days after we had 
finished this paper. C.E.C. is extremely grateful to her for her teaching and en- 
couragement, and mainly for years of friendship. We thank the referee, Dr. an- 
nie Zavagno, for many suggestions that largely improved this presentation. We 
also thank Dr. Y. Yonekura for making his CO data available to us. This project 
was partially financed by the Consejo Nacional de Investigaciones Cientfficas y 
Tecnicas (CONICET) of Argentina under project PIP 5886/05 and PIP 5697/05, 
Agencia PICT 14018, and UNLP under projects 11/G072 and 11/G087. The 
Digitized Sky Survey (DSS) was produced at the Space Telescope Science 
Institute under US Government grant NAGW-2166. 



References 

Brand, J., & Blitz, L. 1993, A&A, 275, 67 

Cappa, C.E., Niemela, VS., Martin, M.C., & McClure-Griffifhs, N.M. 2005, 
A&A, 436, 155 

Carraro, G., Patat, E, & Baumgardt, H. 2001, A&A, 371, 107 
Caswell, J.L., & Haynes, R.F. 1987, A&A, 171, 261 
Claria, JJ. 1977, A&A, 27, 145 

Condon, JJ., Griffith, M.R., & Wright, A.E. 1993, AJ, 106, 1095 

Cesarsky, D., Lequeux, J., Abergel, A., Perault, M., Palazzi, E., Madden, S., & 

Tran, D. 1996, A&A, 315, 305 
Cichowolski, S., Pineault, S., Arnal, E.M., Testori, J.C., Goss, W.M., & Cappa 

C.E. 2001, AJ, 122, 1938 
Comeron, E, Schneider, N., & Russeil, D. 2005, A&A, 433, 955 
Cutri, R.M., Skrutskie, M.F., van Dyk, S., et al. 2003, The IRSA 2MASS All-Sky 

Point Source Catalog, NASA/IPAC Infrared Science Archive 
Deharveng, L., Zavagno, A., & Caplan, J. 2005, A&A, 433, 565 
Deharveng, L., Zavagno, A., Salas, L., Porras, A., Caplan, J., & Cruz-Gonzalez, 

I. 2003, A&A, 399, 1135 
Digel S. W., Hunter S. D., & Mukherjee R., 1995, ApJ, 441, 270 
Dobashi, K., Yonekura, Y, Matsumoto, T. et al. 2001, PASJ, 53, 85 
Duncan, A.R., Stewart, R.T., Haynes, R.E, & Jones, K.L. 1995, MNRAS, 277, 

36 

Dutra, CM., Bica, E., Soares, J., & Barbuy, B. 2003, A&A, 400, 533 
Egan, M. P., Price, S. D., Moshir, M. M., Cohen, M., Tedesco, E., Murdock, T. L., 
Zweil, A., Burdick, S., Bonito, N., Gugliotti, G. M., & Duszlak, J. 1999, The 
Midcourse Space Experiment Point Source Catalog Version 1 .2 Explanatory 
Guide, AFRL-VS-TR- 1999- 1522, Air Force Research Laborartory 
Elmegreen, B.G., & Lada, C.J. 1977, ApJ, 214, 725 
Georgelin, Y. M., Lortet, M. C, & Testor, G. 1986, A&A, 162, 265 
Gum, C.S. 1955, MmRAS, 67, 155 
Hildebrand, R.H. 1983, Q.J.RAS, 24, 267 

Jeffers, H.M., van den Bos, W.H., & Greeby, EM. 1963, Publ. of Lick Obs. 21 
Junkes, N., Fiirst, E., & Reich, W. 1992, A&A, 261, 289 

Lumsden, S.L., Hoare, M.G., Oudmaidjer, R.D., & Richards, D. 2002, MNRAS, 
336, 621 

Martins, E, Schaerer, D., Hillier, D.J., 2005.A&A, 436, 1049 
Mafhys, 1988, A&AS, 76, 427 

McClure-Griffiths, N.M., Dickey, J.M., Gaensler, B.M., & Green, A.J., 

Haverkorn, M., & Strasser, S. 2005, ApJS, 158, 178 
Mezger, P.D., & Henderson, A.P. 1967, ApJ, 147, 471 
Moffat, A.F.J., & Vogt, N. 1975, A&S, 20, 125 

Parker, Q.A., Phillipps, S., Pierce, M.J., et al. 2005, MNRAS, 362, 689 

Prinja, R.K., Barlow, M.J., & Howarfh, I.D. 1990, ApJ, 361, 607 

Rohlfs, K. 1986, "Tools of Radioastronomy", Springer- Verlag Berlin 

Heidelberg, p. 159 
Romero, G.A. 2006, Ph.D. thesis 

Sadler, E.M., & Hunstead, R.W. 2001, ASP Conference Series, 232, 53 
Spitzer, L. 1978, "Physical Processes in the Interstellar Medium", Wiley- 

Interscience publication, New York. 
Vazquez, R.A., & Feinstein, A. 1990, Rev. Mex. AA, 21, 346 
Walborn, N. 1982, AJ, 87, 1300 



12 



C. Cappa, V.S. Niemela, R. Amorin and J. Vasquez,: The environs of Gum 31 



Walborn, N. 1995, in "The r\ Car region: a laboratory of Stellar Evolution", eds. 
Niemela, V., Morrell, N., & Feinstein, A., Rev. Mex. AA, Conf. Series, 2, 
51 

Wood, D.O.S., Churchwell, E. 1989, ApJS, 69, 831 

Yamaguchi, R., Saito, H., Mizuno, N., Mine, Y., Mizuno, A., Ogawa, H., & 

Fukui, Y. 1999, PASJ, 51,791 
Yonekura, Y, Asayama, S., Kimura, K., Ogawa, H., Kanai, Y, Yamaguchi, N., 

Barnes, P.J., & Fukui Y. 2005, ApJ, 634, 476 



